Detecting entanglement of two electron spin qubits with witness operators 
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We propose a scheme for detecting entanglement between two electron spin qubits in a double 
quantum dot using an entanglement witness operator. We first calculate the optimal configuration 
of the two electron spins, defined as the position in the energy level spectrum where, averaged over 
the nuclear spin distribution, 1) the probability to have two separated electrons, and 2) the degree 
of entanglement of the quantum state quantified by the concurrence are both large. Using a density 
matrix approach, we then calculate the evolution of the expectation value of the witness operator 
for the two-spin singlet state, taking into account the effect of decoherence due to quantum charge 
fluctuations modeled as a boson bath. We find that, for large interdot coupling, it is possible to 
obtain a highly entangled and robust ground state. 
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Entanglement - non-classical correlations between 
quantum-mechanical particles - has for a long time been a 
theoretically predicted concept [H without experimental 
proof. The first experimental observation of entangle- 
ment was reported for photon pairs in 1972 and since 
then also entanglement between protons [J kaons [3], 
trapped ions individual neutrons ^ , between an atom 
and a photon [7[ , and between superconducting qubits [8] 
has been demonstrated. An important goal in present- 
day solid-state quantum physics is to generate and de- 
tect (prove) entanglement between individual electrons. 
The motivation behind this quest comes both from the 
fact that entanglement between electrons in a solid-state 
structure has so far not yet been demonstrated and from 
the recent experimental progress in the field of quantum 
information processing in these systems Q, which has, 
among others, led to experimental realization of single- 
and two-qubit manipulations of electron spin qubits in 
quantum dots [lol . 11 1 and coherent control of spins in 



diamond [12 1 



In this Letter we propose a scheme for detecting entan- 
glement in the former system, namely entanglement of 
two electron spin qubits in a double quantum dot. Many 
aspects of this quantum system, such as hyperfine cou- 
pling to the nuclear spins (IS-lSj , spin blockade ,l(i JJ'] , 
and effects of applying a slanting magnetic field [18| are 
currently active topics of research. Our proposal con- 
sists of preparing and tuning the two-spin system such 
that the ground state contains a large entangled com- 
ponent and measuring a so-called entanglement witness 
operator [l9| to demonstrate the presence of this entan- 
glement. Entanglement witnesses are hermitian opera- 
tors that are designed to detect a specific entangled state. 
Their expectation value is positive for all separable quan- 
tum states (in the class of states considered) and negative 
for at least one entangled state, usually the state the ex- 
periment aims to create. Entanglement witnesses have 
been used to detect entangled states in trapped ions sys- 
tems or entangled photons states [20l |. 



Theoretical proposals for detecting entanglement be- 
tween individual electrons using witness operators are 
scarce. Ref. [2l| presents a proposal for implementing 
witness operators to detect electron-hole entanglement 
in multiterminal conductors in the presence of noise due 
to random accumulated phases. In a previous work 22], 
we have proposed a turnstile mechanism as a suitable set- 
up for demonstrating entanglement between two electron 
spins in a double quantum dot, assuming phenomenolog- 
ical decoherence times for spin relaxation (Ti-time) and 
decoherence (T2-time). In the present paper we focus 
on the dynamic evolution of the entanglement witness 
operator under the influence of decoherence due to quan- 
tum charge fluctuations, which are caused by changing 
gate voltages to control the system. Starting from the 
two-electron Hamiltonian, and assuming a large exter- 
nal magnetic held, we first calculate the energy levels 
and eigenstates in the three-dimensional Hilbert space 
spanned by the |S'(0,2)), |S'(1, 1)) and |To(l, 1)) eigen- 
states, where |<S'(n, m)) (|To(n, m))) represents the singlet 
(triplet) state with n electrons in the left and m electrons 
in the right dot. We then average over the nuclear field 
components and calculate the average probability of hav- 
ing one electron in each dot (Pii)nuci and the concurrence 
(C')nuci of the entangled component of the ground state. 
This allows us to determine the "optimal" configuration 
in the level diagram in which both of the latter quantities 
are large. We then investigate the effects of decoherence 
of the entangled quantum state due to quantum charge 
fluctuations, modeled as a boson bath. Using the Born- 
Markov approximation, we calculate the reduced density 
matrix p[t) of the system and the evolution of the expec- 
tation value of the witness operator Tt:{W p{t)). We flnd 
the optimal values in our parameter space that provide a 
highly entangled ground state which is also robust under 
the action of charge fluctuations. 

Model. - We consider a charge configuration of the 
double quantum dot where each electron can be located 
in a different quantum dot (1,1), or both of them in the 
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right dot (0,2). In the (1,1) configuration the four acces- 
sible states are the singlet 15(1, 1)) and the three triplet 
states \Ti{l, 1)). In the (0,2) configuration only the sin- 
glet |S'(0, 2)) can be populated, the triplet states \Ti{0, 2)) 
having much higher energies. The spin-preserving part 
of the Hamiltonian is given by: 



Ho = A|5o2)(S'o2| +t(|S'o2)(5ii| + |5ii)(S'o2|). (1) 

The tunneling parameter t couples the two singlets allow- 
ing one electron to transfer between the two dots, while A 
is the energy difference between the (0,2) - and the (1,1) - 
singlets. Both parameters can be externally controlled by 



changing gate voltages and are used to tune the system 
to the desired configuration. In addition, the hyperfine 
interaction of the electrons with the nuclear spins mixes 
the singlet and triplet states in the (1,1) configuration. 
Each electron interacts with a large number of nuclear 
spins in the left (L) and right (R) dot. The global action 
of these nuclear spins can be included in a single operator 

R which allows us to treat the hyperfine interaction as 
the interaction between the electrons and this apparent 
magnetic field. We assume that this nuclear magnetic 
field remains unchanged over the typical timescale of the 
electron spin evolution 25] and hence it can be treated 
as a classical magnetic field. The combined action of the 
nuclear magnetic field and of an external magnetic field 
Bext is given by 



I 

= B:{\T+){T+\ - |T_)(r_|) + [B:\To){S,^\ + ^l±^\To){T±\ + ±^l^^\S^,){T±\ + H.c.) , (2) 



where = (B^ - B^)/2, B, = (B^ + B^)/2 + B.^t z. 
The total Hamiltonian of the system is then given by 
Hdqd ~ Hq + Hspin 17], and its dynamics can be ex- 
ternally controlled by tuning A, t and B^xt- Assuming a 
large external magnetic field B^xt the triplet states \Ti) 
are split off by the Zeeman energy, and the state space 
is reduced to {\To), \S(1, 1)), 15(0,2))}. Under these cir- 
cumstances it is possible to derive an effective Hamilto- 
nian lisl. 



H,ff = A\So2){So2\+t\So2){Sn\+M\Sn){To\+li.c. (3) 

where M = B^ + [{B^ + iBy){B^ - iBy) + c.c] /2BI. 
The eigenstates of the Hamiltonian ([3]) are given by 



M 

Ei 



\To) + \Sn) + E^\So2) 



(4) 



(B.-A)^ 



where Ei are their corresponding energies {Eq < Ei < 
i?2)- These energies are given by i^oCi) = A/3 — (si + 
S2)/2 + (-)i\/3(£]_ - S2)/2 and E2 = A/3 + si + S2, with 
si,2 = (r± v/g3 + r2)i/3^ r = A(-18M2 -h Qt^ 2A2)/54 
and q = -(3(M2 -h t^) + A2)/9 and are plotted in 

Fig.m 

The entangled component of the ground state is given 
by its projection onto the (l,l)-subspace, 



,l^) = {l + MyE',)-'^'(^fjTo) + \Sn)y (5) 

The amount of entanglement of |V'o^) can be quantified 
by the concurrence 2j|. Both the concurrence C and 



the probability Pn of having one electron in each dot 
must be large to have a highly entangled state. Their 
mathematical expressions are given by 
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(6) 



(£o-A) 



The concurrence C is maximal when |i?o| ^ M, while 
Pii is maximal in the presence of a large detuning (A 3> 
|£'o|). It can be observed in Fig. [T]that both conditions 
can only be fulfilled for intermediate values of A where 
A > |£'o| > M (typically M - 0.1- 1 ^leV). To find the 
optimal A and t values we consider the maximization of 
their product ^ = PuC. 

Averaging over the nuclear field distribution. - Given 
the random character of the nuclear magnetic field, 
any observable can be calculated taking the average 
over the probability distribution that characterizes B^^. 
This distribution is Gaussian with variance (B^^) = 
El/Neff, with En « 0.135 meV for GaAs and Neff « 
10^ for typical dots [13, [l^- For external magnetic fields 
Bext > 2.5 ^eV (100 mT) the transversal components 
of the nuclear magnetic field (B^y) are negligible and 
the hyperfine interaction parameter M is mainly given 
by the difference of the z components of the nuclear 
magnetic fields (Af ~ B^). The probability distribu- 
tion for M is then approximately Gaussian with variance 
a' = (B^^^,)/2. 

Fig. [2] shows the averaged success rate {£,)nuci ob- 
tained using the previous approximated Gaussian distri- 
bution for M. The entanglement of the ground state 
is low for A < where Pu < 0.5, and for small val- 
ues of the interdot coupling t where the ground state 
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FIG. 1. (Color online) Energy level diagram of the double dot 
system for interdot coupling t = 5 fieV. 



-150 150 300 450 600 750 

A (/leV) 



is nearly a separable state. The black line links the 
optimal values {Aopt,topt} that maximize {C)nuci- For 
topt > 5 ^ieV, and in a large region nearby, the ground 
state is highly entangled {{£,)nuci > 0.95) up to values of 
0.99 for topt = 20 jieV. Obtaining the average over the 
nuclear probability distribution is usually a difficult task, 
but in this case the value of the hyperfine parameter M in 
Eqs. ^ can to a good approximation be substituted by 
the standard deviation of its Gaussian probability dis- 
tribution mM))nuci ~ ^(vW)) = ^(a)). We have 
checked that this approximation holds in the highly en- 
tangled region and only fails for small t and large A. In 
that region i?o — M, and small variations in M have a 
great impact on the (small) value of the concurrence C, 
and hence on ^. 

The entanglement of the system can be detected us- 
ing witness operators. For an entangled state \ip) — 
7i|To) + 72|S'(1, 1)) (with 72 > 71) the optimal operator 
to detect its entanglement is given by W = \Tq) (TqI -|- 
\S{1, 1)) (5(1, 1)1 + |To) (5(1, 1)1 + |5(1, 1)) (Tol 0. Its 
expectation value is given hy Tr{W\ip){(p\) = (71— 72)/2, 
which is proportional to the concurrence {Tr{W\ip) {ip\) = 
—C/2). This witness can be rewritten as W = 1/2 — 
|5(1, 1))(5(1, 1)|. In order to measure the expectation 
value of VF, it is then enough to measure the probability 
P[5(l,l)] of the ground state to be the |5(1,1)) state. 
In a double quantum dot this probability can easily be 
measured using a quantum point contact (QPC) [9j. Our 
proposal to detect the entanglement is then the following: 
To initialize the system in the |5(0, 2)) state (point A in 
Fig. [2]) and adiabatically increase A to the optimal po- 
sition {Aopt,topt} (point B in Fig. [2]). A QPC can then 
be used to measure P(S(1,1)) and obtain the expectation 
value of the witness operator. 

Time evolution of the witness operator. - When mod- 
ifying gate voltages to drive the system to its optimal 
position and to measure the probability P[5(l,l)], the 
resulting charge fluctuations cause the tunnel coupling t 
and energy offset A to fluctuate. The inset of Fig. [3] 

lOmK and topt > 5^eV, 



FIG. 2. (Color online) Averaged entanglement probability 
< C >nuci- The black line represents the optimal configura- 
tion {Aopt,topt} that maximizes the entanglement of the sys- 
tem. The white solid line represents the trajectory followed 
from the initial state A to the final state B. 



(i.e. in the highly entangled region) so that quantum 
fluctuations are important [2^. In the last part of this 
Letter, we explore the influence of these fluctuations on 
the time evolution of the expectation value of the entan- 
glement witness. 

Modelling the environment as a bosonic bath [27|, the 



Hamiltonian of the double quantum dot plus the envi- 
ronment is given by 



H = Hdqd + VtAt + VaAa + Hi 



bath.t 



H, 



bath, A 1 



(7) 



shows that Egap > ksT for T 



where = |5(1, 1))(5(0, 2)| + |5(0, 2))(5(1, 1)|, Va = 
|5(0,2))(5(0,2)|, = T.k<^^Ablk + Kk), e^ndHbath,^ = 

^kfvjJk{b\^. + bi^k), with i — t, A. The bosonic 
baths are characterized by symmetric and antisym- 
metric spectral functions 5^(w) which are related by 
S^{uj) — coth{huj/2kbT)S~ (uj). We assume a bath with 
Lorentzian damping S~{uj) = ct^^^j^^jjj^JY ' '^^^^ a 
high cut-off frequency. We also assume weak coupling 
between the system and the bath, and short bath cor- 
relation times. Following the Bloch-Redfield approxima- 
tion, the time-dependent reduced density matrix, writ- 
ten in the basis of the eigenstates of Hdqd is given by 

dtPab = -iUJabPab + J2cd RabcdPcd, whcrC tvjJab = Ea - Eb 

and Rabcd is the Bloch-Redfield tensor [27|. The sum 
in this equation extends over terms with uJab — ^cd <C 
1/Ar (the so-called secular constraint), where At is the 
timescale of the Markovian course-grained evolution [2I] . 
Considering the initial condition p(0) = |-!/'o)(V'o| the 
only relevant components of the Bloch-Redfield tensor 
are i?oooo = —^1100 and i?iiii = — i?ooii- The time evo- 
lution for the populations of the ground and first excited 
states are then given by 
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FIG. 3. (Color online) Time evolution of the expectation 
value of the witness operator for T = 10 mK, a = 0.007 [2^ 
and topt ~ 20 fieV (orange line), 10 iJ,eV (green line), 5 f^eV 
(black line). Dashed lines mark their asymptotic values. In- 
set: Energy gap between ground and first excited states in 
the optimal configuration {Aopt,topt} (blue line), the red line 
indicates the thermal energy fcsT — 0.86 fieV . 



Poo(t) 
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(8) 
(9) 



Rewriting the eigenstates given in ^ as \ipi) — ailTo) + 
PilSii) + 7i|S'o2), the coefficients Rujj for the optimal 
values {Aopt,topt} are approximately given by Rnjj ~ 
2G{ujji){l3aj + Pjlif/^^ with G(a;„) = 2T:^ah^uj[l + 
coth(fta;/2fcsr)]. Introducing these expressions in Eq. 
(|8]) and ^ we can obtain the time dependent probabil- 
ity P[5(l, l)](r) = 1(5(1, l)|p(r)|5(l, 1))| which provides 
the expectation value of the entanglement witness plot- 
ted in Fig. 131 For these optimal values P[S'(1,1)] is 
approximately given by 



P[5(l,l)] «Poo(r)«i/3o72. 



(10) 



As can be seen in Fig. [3l the larger the value of the in- 
terdot coupling t the longer the system retains a large 
amount of its initial entanglement. If the ground state 
is in its optimal configuration {lS.opt,topt\ the loss of en- 
tanglement caused by the quantum charge fluctuations in 
the environment becomes important for timescales of the 
order of microseconds |25|] . These times are much longer 
than the timescales required to manipulate electrons con- 
fined in a double quantum dot, typically of the order of 
1-100 ns We note that for a given topt the decay of 
the entanglement witness expectation value is smaller for 
A < Aopj. The loss of entanglement of the system while 
it is adiabatically moved to its optimal configuration is 
then even smaller that what is seen in Fig. |3l 

Conclusion. - We have proposed a scheme to detect the 
entanglement between two electrons in a double quantum 



dot and shown that it is possible to maximize the entan- 
glement between the two electrons in the ground state 
of the system. This entanglement can be detected using 
an entanglement witness, and its expectation value can 
be easily measured using a quantum point contact. We 
have also found that both in the optimal configuration 
that maximizes the entanglement and in its trajectory to 
this configuration, the entanglement of the ground state 
is robust against charge fluctuations in the environment. 

We would like to thank Fabian Bodoky for useful dis- 
cussions. This work was supported by the EU project 
MIDAS. 
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